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Abstract 

Background: The incorporation of histone variants into nucleosomes is one of the main strategies that the cell 
uses to regulate the structure and function of chromatin. Histone H2A.Z is an evolutionarily conserved histone H2A 
variant that is preferentially localized within nucleosomes at the transcriptional start site (TSS). H2A.Z reorganizes the 
local chromatin structure and recruits the transcriptional machinery for gene activation. High expression of H2A.Z 
has been reported in several types of cancers and is causally linked to genomic instability and tumorigenesis. 
However, it is not entirely clear how H2A.Z overexpression in cancer cells establishes aberrant chromatin states 
and promotes gene expression. 

Results: Through integration of genome-wide H2A.Z ChlP-seq data with microarray data, we demonstrate that 
H2A.Z is enriched around the TSS of cell cycle regulatory genes in bladder cancer cells, and this enrichment is 
correlated with the elevated expression of cancer-promoting genes. RNAi-mediated knockdown of H2A.Z in the 
cancer cells causes transcriptional suppression of multiple cell cycle regulatory genes with a distinct decrease in 
cell proliferation. H2A.Z nucleosomes around the TSS have higher levels of H3K4me2/me3, which coincides with 
the recruitment of two chromatin factors, WDR5 and BPTF. The observed recruitment is functional, as the active 
states of H2A.Z target genes are largely erased by suppressing the expression of WDR5 or BPTF, effects resembling 
H2A.Z knockdown. 

Conclusions: We conclude that H2A.Z is overexpressed in bladder cancer cells and contributes to cancer-related 
transcription pathways. We also provide evidence in support of the engagement of H3K4me2/me3 and WDR5/BPTF 
in H2A.Z-induced cancer pathogenesis. Further studies are warranted to understand how H2A.Z overexpression 
contributes to the recruitment of the full repertoire of transcription machinery to target genes in bladder cancer cells. 

Keywords: Histone, H2A.Z, Chromatin, Bladder cancer, ChlP-seq, Microarray, Gene expression 



Background 

Gene expression in eukaryotes takes place in the context 
of chromatin, which essentially acts as a physical barrier 
to the binding of transcription factors to certain genomic 
regions. The basic unit of chromatin is the nucleosome, 
which consists of approximately 147 bp of DNA wrapped 
around an octamer of four core histones (H2A, H2B, 
H3, and H4) [1]. Despite its high degree of compaction, 
chromatin structure is dynamically modulated to allow 
access of genomic DNA during vital cellular processes. 
Several remodeling processes that alter chromatin 
structure have been identified and categorized based on 
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their modes of action. The three most prominent processes 
are histone modification, histone variant exchange and 
adenosine triphosphate (ATP)-dependent chromatin 
remodeling [2-5]. Histone variants are nonallelic isoforms 
of canonical core histones that have been identified 
primarily from the H2A and H3 families in numerous 
organisms. Although expressed in smaller amounts 
than canonical histones, histone variants are involved 
in a variety of cellular activities with specific genomic 
localization patterns. 

H2A.Z is an essential variant of the histone H2A and 
has distinct functions in regulating chromatin dynamics. 
At the amino acid sequence level, H2A.Z shares about 
60% identity with canonical H2A in mammalian cells 
[6]. As for other histone variants, the replacement of 
histone H2A with H2A.Z is replication-independent and 
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is catalyzed by ATP-dependent exchange factors. H2A.Z 
preferentially localizes at the transcription start site 
(TSS) where it frequently flanks nucleosome-deficient 
regions [7-11]. The presence of H2A.Z at the -1 and +1 
nucleosomes adjacent to the TSS has been linked to 
dynamic changes in gene expression [12-15], but the 
underlying mechanism is not fully understood. Since 
H2A.Z-containing nucleosomes wrap DNA more weakly 
than canonical H2A nucleosomes, they are more suscep- 
tible to nuclease digestion and salt-induced destabilization 
[16,17]. This structural instability is most likely driven 
by amino acid substitutions at the interface between 
H2A.Z and H3/H4 [18]. Consistent with its effects on 
gene transcription, the presence of H2A.Z correlates with 
active histone modifications and RNA Pol II occupancy at 
gene promoters, which in turn influence transcriptional 
output [7,13,19]. These observations are in agreement with 
the concept that H2A.Z creates an active transcription 
environment by altering DNA accessibility in chromatin 
and facilitating the recruitment of transcription regulators. 
Of special relevance to the present study, overexpression 
and redistribution of H2A.Z have been noted as epigenetic 
processes to stimulate the development and progression 
of certain types of cancer [20-23] . The oncogenic function 
of H2A.Z is mediated through increased cell proliferation 
caused by cell cycle deregulation. Therefore, the important 
questions to address are how H2A.Z maintains chromatin 
in a certain state, and whether H2A.Z overexpression 
observed in cancer cells drives the inappropriate activation 
of cancer-related genes. 

In this work, we exploited microarray and chromatin 
immunoprecipitation sequencing (ChlP-seq) to examine 
the impact of H2A.Z overexpression on transcriptional 
states of bladder cancer cells, and identified groups of 
genes, the dysregulation of which is linked to H2A.Z 
enrichment around the TSS. We also found that H2A.Z 
nucleosomes are enriched for the active histone modifica- 
tion H3K4me2/me3 that is likely responsible for recruiting 
WDR5 and BPTF to H2A.Z target genes. 

Results 

H2A.Z is overexpressed in bladder cancer cells and 
stimulates cell proliferation 

As H2A.Z expression is often misregulated in cancer 
cells, we first examined H2A.Z levels in human bladder 
cell lines by Western blot. H2A.Z was expressed at much 
higher levels in the six bladder cancer cell lines (SCaBER, 
J82, LD611, RT4, HT1376 and T24) compared to the 
normal bladder cell line UROtsa (Figure 1A). Another 
demonstration of the relationship between H2A.Z over- 
expression and tumorigenesis came from the results 
obtained from immunohistochemical analysis of tissue 
microarrays containing 16 bladder cancer samples and 
8 normal control samples. As illustrated in Figures IB 



and 1C, the staining for H2A.Z was strong or moderate 
in the majority of bladder tumor tissues (approximately 
90%), while the staining was weak or negative in most 
normal tissues (approximately 85%). 

We next checked whether the high-level expression of 
H2A.Z in cancer cells contributes to cell growth over a 
period of four days. The 3-4, 5- (dimethyl- thyazol-2-yl) -2, 
5-diphenyltetrazolium (MTT) assays reproducibly showed 
that LD611 cancer cells expressing high levels of H2A.Z 
grow much faster compared to normal UROtsa cells 
(Figure ID). These findings were further corroborated 
by the results showing that the ectopic expression of 
H2A.Z led to a significant increase in the growth of 
UROtsa cells (Figure IE and Additional file 1: Figure 
SI). To evaluate the effects of H2A.Z overexpression 
more directly, endogenous H2A.Z was depleted in LD611 
cells by using a lentiviral shRNA infection system, and 
changes in cell growth were analyzed. Western blotting 
confirmed that the stable transfection of H2A.Z 
shRNA efficiently silenced H2A.Z expression in the 
cell (Additional file 1: Figure S2). The MTT assays 
clearly showed that suppressing H2A.Z expression 
decreased the growth of LD611 cells (Figure IF). 
These observations are consistent with the hypothesis that 
H2A.Z is one of the key players governing proliferation 
of bladder cancer cells. 

H2A.Z shows a unique localization pattern in bladder 
cancer cells 

To gain mechanistic insights into the observed effects of 
H2A.Z, we conducted ChlP-seq of H2A.Z in the LD611 
cancer cell line. The UROtsa cell line was used as a normal 
control to identify the genomic redistribution of H2A.Z 
in the cancer cell line. We compared the average ChIP 
enrichment signals over the regions spanning 10 kb 
around the TSS of all human genes as detailed in the 
Methods. As expected, H2A.Z occupied nucleosomes 
surrounding the TSS in both normal and cancer cells 
(Figure 2A). However, the enrichment of H2A.Z was 
more pronounced in the bladder cancer cells than in the 
normal cells. Heatmap analysis with k-means clustering of 
ChlP-seq datasets also indicated that H2A.Z nucleosomes 
were more highly enriched in the immediate upstream 
and downstream of the TSS in the cancer cells (Figure 2B). 
To distinguish the H2A.Z enrichment from the genomic 
background, the expectation maximization-based curve 
fitting method that uses the normalized locus-specific 
chromatin state (NLCS) was employed as detailed previ- 
ously [24]. Based on the log 2 transformed NLCS values 
and false discovery rate (FDR), we found that 750 and 200 
genes were enriched for H2A.Z in the cancer and normal 
cells, respectively. To identify the genes involved in 
bladder cancer development, we further classified the 
cancer-specific genes into the four distinct groups 
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Figure 1 H2A.Z overexpression and bladder cancer cell growth. (A) Normal (UROtsa) and bladder cancer (SCaBER, J82, LD61 1, RT4, HT1376 
and T24) cells were lysed in RIPA buffer, and equal amounts of lysates were analyzed by Western blot using anti-H2A.Z antibody. Actin served as 
a control for equal protein loading. (B) Tissue microarrays containing 16 bladder tumor and 8 normal tissue samples were subjected to 
immunohistochemistry with H2A.Z antibody. Representative high-powered images are shown. (C) Immunostaining scores of H2A.Z in normal 
and bladder cancer tissues. The graph indicates the percentage of sections with different scores (negative, weak, moderate and strong). (D) 
Proliferation of UROtsa and LD61 1 cells was determined by MTT colorimetric assays at the indicated time points. Data shown represent mean ± SD 
from three replicates. (E) UROtsa cells were transfected with H2A.Z as in Additional file 1 : Figure S1, and cell proliferation was measured by MTT assay 
over a period of four days. All reactions were performed in triplicate. (F) LD61 1 cells were mock-depleted or depleted of H2A.Z as in Additional file 1 : 
Figure S2, and cell proliferation was monitored by fvuT assays over a period of four days. IVHT assay, 3-4, 5-(dimethyl-thyazol-2-yl)-2, 
5-diphenyltetrazolium assay. 



based on their FDR values and NLCS. A total of 480 
out of 750 genes were classified into the groups 1, 2 and 3 
with P <0.05, suggesting that these genes are specifically 
enriched for H2A.Z in the bladder cancer cells (Figure 2C 
and Additional file 2: Table SI). To validate ChlP-seq 
results, we chose three H2A.Z-enriched genes in the 
normal and cancer cells and carried out conventional 
ChIP assays. Expectedly, H2A.Z was detected around 
the TSS of these genes in both cell lines (Additional file 1: 
Figure S4). Ingenuity pathway analysis also identified a set 
of genes related to bladder cancer signaling as well as 
DNA repair (Figure 2D). 

H2A.Z redistribution correlates with increased expression 
of proliferation-related genes 

The aberrant enrichment of H2A.Z within nucleosomes 
adjacent to the TSS in the bladder cancer cells suggests 
that H2A.Z may establish distinct transcription states. 
To interrogate this possibility, we conducted gene expres- 
sion profiling with total RNA isolated from the normal 
and cancer cell lines. Our detailed comparison of the 



transcriptional profiles identified a large number of differ- 
entially expressed genes in the bladder cancer cells, with 
2,115 genes being upregulated and 2,070 genes being 
downregulated (Figure 3A and Additional file 3: Table 
S2). Quantitative reverse transcription polymerase 
chain reaction (qRT-PCR) analyses of six upregulated, 
six downregulated and six unaffected genes validated 
our microarray results (Additional file 1: Figure S5). 
To infer which of the genes showing altered expres- 
sion in the cancer cells were potentially regulated by 
H2A.Z, we combined the ChlP-seq and gene expres- 
sion data. Among the 2,115 genes upregulated in the 
cancer cells, 133 genes exhibited higher levels of H2A.Z. 
Conversely, only 23 genes out of the 2,069 downregulated 
genes showed H2A.Z enrichment (Figure 3A and Add- 
itional file 4: Table S3). In light of the fact that the ac- 
quisition of H2A.Z around the TSS is mainly correlated 
with gene activation, we decided to focus on the 
upregulated genes. Functional annotation of differen- 
tially expressed genes using ingenuity pathway analysis in- 
dicated that genes involved in cell growth and proliferation 
are significantly enriched in this category (Figure 3B). 
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Figure 2 H2A.Z localization patterns in bladder cancer cells. (A) Average chromatin immunoprecipitation (ChIP) enrichment signals are 
shown over regions spanning 10 kb around the transcription start sites (JSSs) of all the human genes from UCSC Refseq database. Green and red 
lines indicate the level of H2A.Z enrichment in normal UROtsa and cancer LD61 1 cells, respectively. (B) Tag density plots display H2A.Z 
distribution relative to the TSS. Several representative peaks specific for LD61 1 and UROtsa cells were also shown. Each horizontal line represents 
an individual gene, and the enrichment values for H2A.Z are illustrated by color density. (C) H2A.Z-enriched genes in LD61 1 cells were classified 
into four distinct groups based on false discovery rate (FDR) and the level of locus-specific chromatin state (NLCS) density. Groups 1 and 2 show 
NLCS ratio >8.5 and 8 to 8.5 with FDR <0.01 and FDR <0.05, respectively. Group 3 has FDR <0.05 with NLCS level from 7.5 to 8.5. Group 4 has 
FDR >0.05 with NLCS <8. (D) H2A.Z-enriched genes in LD61 1 cells were subjected to ingenuity pathway analysis. Seven top biological functions 
including bladder cancer signaling were identified with a threshold set at P = 0.05. 



When we analyzed a group of the positively regulated 
proto-oncogenes (BIRC3, CCND, DDX43, FADD, MEN1, 
NTSE, POLDIP2, PPFIA1, STK3, and YAP1), we detected 
a robust signal for H2A.Z around the TSS (Figure 3C, 
LD611). In the normal bladder cells, these genes 
displayed much lower levels of H2A.Z (UROtsa). To 
further validate the ChlP-seq and gene expression 
data, we selected the 10 active genes from the H2A.Z 
regulon and performed conventional ChIP assays. 
Expectedly, H2A.Z enrichment was evident in these 
genes in the bladder cancer cells, but similar assays in 
the normal cells showed minimal accumulation of 
H2A.Z (Figure 3D). Moreover, there was a significant 
correlation between H2A.Z enrichment and gene acti- 
vation, as evidenced by qRT-PCR analysis (Figure 3E). 
To the best of our knowledge, this is the first report 



linking H2A.Z enrichment to active transcription 
states in bladder cancer cells. 

WDR5 and BPTF preferentially interact with H2A.Z 
nucleosomes 

To study how H2A.Z directs oncogenic transcription in the 
bladder cancer cells, we decided to identify proteins that 
selectively interact with H2A.Z nucleosomes. Accordingly, 
the LD611 cells were transfected with plasmids expressing 
Flag/HA-tagged versions of H2A and H2A.Z. After con- 
firming the comparable expression of ectopic H2A.Z and 
H2A (data not shown), soluble chromatin was prepared 
from the transfected cells and digested with micrococcal 
nuclease to yield mainly mononucleosomes. Mononuc- 
leosomes containing ectopic histones and their bound 
proteins were purified by sequential immunoprecipitations 
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Figure 3 Effects of H2A.Z enrichment on gene expression. (A) Venn diagram shows H2A.Z-enriched genes that are upregulated (133 genes) 
and downregulated (23 genes) more than 1.5-fold with a P value of <0.05 in LD61 1 bladder cells. (B) Functional groups associated with the 133 
genes that are enriched by H2A.Z and are upregulated in LD61 1 cancer cells. Biological processes are graphically depicted by P value. (C) 
Representative Integrative Genomics Viewer (IGV) genome browser tracks show the levels of H2A.Z near the TSS in LD61 1 (red) and UROtsa (blue) 
cells. Maximum sequencing depth over the displayed area is indicated on the right side in parentheses along with name of the genes. Different 
genomic coordinates and genome window size for BIRC3 (chrl 1:101,691,391-101,717,344; 26 kb); CCND1 (chrl 1:69,163,054-69,180,423; 17 kb); 
DDX43 (chr6:74,1 59,006-74,1 86,010; 26 kb); FADD (chrl 1 :69,724,91 7-69,733,1 56; 8 kb); MEN1 (chrl 1 :64,325,562-64,337,342; 1 1 kb); NT5E 
(chr6:86,21 4,021 -86,264,228; 50 kb); POLDIP2 (chrl 7:23,695,786-23,71 0,730; 15 kb); PPFIA1 (chrl 1 :69,792,454-69,91 0,255; 118 kb); STK3 
(chr8:99,505,659-99,937,462; 432 kb) and YAP1 (chrl 1:101,484,402-101,61 1,364; 127 kb) are shown along with Reference Sequence (RefSeq) 
data. (D) LD61 1 and UROtsa cells were analyzed by quantitative chromatin immunoprecipitation (ChIP) using H2A.Z antibody to validate 
ChlP-sequencing results. Immunoprecipitated DNA was amplified by quantitative PCR (qPCR) with primers specific for the 10 genes. The primers 
are listed in Additional file 6: Table S5. Percent input was determined as the amount of immunoprecipitated DNA relative to input DNA. Error 
bars represent standard deviation of three independent experiments. (E) The RNA samples from LD61 1 and UROtsa cells were analyzed by 
quantitative reverse transcription PCR (qRT-PCR) using primers listed in Additional file 6: Table S5. Expression level was normalized to that of actin, 
and average and standard deviation of three independent experiments are shown. 



with anti-Flag and anti-HA antibodies under stringent 
conditions (300 mM NaCl and 0.1% NP40). Coomassie 
brilliant blue (CBB) staining confirmed the stoichiometry 
of H2B, H3, H4 and ectopic H2A/H2A.Z in the purified 
nucleosomes (Figure 4 A, CBB). Western blotting with 
H2A and H2A.Z antibodies also indicated that ectopic 
H2A and H2A.Z were incorporated into nucleosomes 
at levels substantially higher than their endogenous 
counterparts (a-H2A and a-H2A.Z). Silver staining of 



an SDS-polyacrylamide gel containing the proteins that 
were co-purified with H2A or H2A.Z nucleosomes revealed 
the presence of multiple extra bands specific to H2A.Z 
nucleosomes. To identify these extra polypeptides, the puri- 
fied proteins were subjected to liquid chromatography- 
tandem mass spectrometry (LC-MS/MS). The data 
were filtered by excluding any proteins observed in 
control pull downs with Flag-HA-tag. This approach 
identified 27 proteins specific to H2A.Z nucleosome 
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Figure 4 Preferential binding of WDR5 and BPTF to H2A.Z nucleosomes. (A) LD61 1 cells were transfected with H2A and H2A.Z expression 
vectors, and mononucleosomes were prepared by MNase digestion as recently described [25]. Mononucleosomes containing ectopic H2A 
(lane 1) and H2A.Z (lane 2) were sequentially immunoprecipitated from total mononucleosomes with anti-Flag and anti-HA antibodies. Histone 
compositions of the purified nucleosomes were analyzed by 15% SDS-PAGE followed by Coomassie brilliant blue (CBB) staining. The relative 
levels of endogenous/ectopic H2A and H2A.Z in the purified nucleosomes were determined by Western blotting using H2A and H2A.Z antibodies. 
(B) The proteins co-purified with H2A and H2A.Z nucleosomes were resolved in 4 to 20% SDS-polyacrylamide gels, and then visualized by silver 
staining. (C) Venn diagram depicts the number of proteins associated with H2A nucleosomes (65), H2A.Z nucleosomes (85) and those in common (58). 
(D) The proteins specifically co-purified with H2A.Z nucleosomes were functionally classified by gene ontology analysis. The number of proteins in 
each functional group is shown. (E) The preferential binding of WDR5 and BPTF to H2A.Z nucleosomes were analyzed by Western blotting. (F) 
Mononucleosomes containing ectopic H2A or H2A.Z were purified as in (A), and subjected to Western blotting with antibodies specific for the 
indicated histone modifications. 



pull downs (Figure 4C and Additional file 5: Table 
S4). Gene ontology analysis of the proteins revealed 
a striking overrepresentation of proteins related to 
chromatin reorganization and transcriptional control 
(Figure 4D). 

A noteworthy observation emerged from our purifica- 
tion was substantial interactions of H2A.Z nucleosomes 
with WDR5 and BPTF, which are known to interact with 
H3K4me tails and play important roles in regulating 
chromatin function [26-28]. The observed interactions 
are specific, as H2A nucleosomes failed to interact with 
these two factors (Figure 4E). Further analyses revealed 
that the overall levels of neither H3K9me3 nor 
H3K27me3 were altered in H2A.Z nucleosomes (Figure 4F, 



a-H3K9me3 and a-H3K27me3). However, in the same 
nucleosome samples, H2A.Z nucleosomes showed much 
higher levels of H3K4me2/me3 as compared to H2A nu- 
cleosomes (a-H3K4me2 and a-H3K4me3). These results 
underscore the notion that the observed binding of 
WDR5 and BPTF depends on their ability to recognize 
H2A.Z-facilitated H3K4me2/me3 in the nucleosome. 

H2A.Z-mediated gene activation requires WDR5 and BPTF 

The fact that WDR5 and BPTF co-purified with H2A.Z 
nucleosomes raised the question of whether they play 
a potential role in regulating H2A.Z target genes. To 
address this question, we individually knocked down 
H2A.Z, WDR5, and BPTF in LD611 cancer cells and 
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looked at the expression of H2A.Z target genes (Additional 
file 1: Figure S2). As expected, qRT-PCR analysis showed 
that knockdown of endogenous H2A.Z inactivated the 
10 representative target genes (Figure 5A). These genes 
were also expressed at much lower levels in cells de- 
pleted of either WDR5 or BPTF as compared to con- 
trol cells, indicating that WDR5 and BPTF are 



critical for transcriptional activation of the target 
genes. The requirements for WDR5 and BPTF were 
highly specific, as the two control genes, CASP8 and 
SRCAP, were minimally affected by silencing endogen- 
ous WDR5 and BPTF. 

Because WDR5 and BPTF are necessary for H2A.Z to 
upregulate gene transcription in the bladder cancer cells, 




B 



CD 
O 

A 
>— 
O 
(/) 
.Q 
< 



NCsh 



H2AZsh WDR5sh 



BPTFsh 



1.51 



1.0" 



0.5 • 



NCsh 
H2A.Zsh 
WDR5sh 
BPTFsh 





Bar, 50 u.m 



T- 

0 



T 1 

1 2 

Time (day) 



T- 

3 



4 



£ 120- 

| 100- 

| 80- 

I 60- 

3 40- 

I 20- 

DC o-l 



I 

lla 



/ y / f 



4 



Figure 5 Cooperative function of H2A.Z 7 WDR5 and BPTF. (A) After silencing H2A.Z, WDR5 or BPTF in LD61 1 cells, the mRNA levels of the 10 
genes, which show higher H2A.Z levels around the transcription start site (TSS) and augmented expression, were quantified by qRT-PCR. The 
values are expressed as fold changes from the levels in undepleted cells. (B) IVHT proliferation assays were carried out in triplicate using LD61 1 
cells depleted of H2A.Z, WDR5 or BPTF as indicated. (C) LD61 1 cells depleted of H2A.Z, WDR5 or BPTF were seeded onto the upper chamber, and 
then allowed to migrate toward 10% fetal bovine serum (FBS) in the lower chamber. The migrated cells were photographed and counted. The Y 
axis indicates the relative percentage of the migrated cell. Three independent experiments in triplicate wells were performed. IvuT assay, 3-4, 
5-(dimethyl-thyazol-2-yl)-2, 5-diphenyltetrazolium assay. 
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we also check whether they are functionally linked to 
H2A.Z-induced cell growth. As summarized in Figure 5B, 
when H2A.Z was depleted in the LD611 cancer cells, the 
growth of the cancer cells was significantly inhibited 
over a four-day period. More importantly, knockdown 
of WDR5 or BPTF had the same inhibitory effect as 
knockdown of H2A.Z. In accordance with these results, 
individual knockdown of H2A.Z, WDR5 and BPTF also 
reduced cell migration severely, displaying about 70% 
lower migration than control cells (Figure 5C). Considered 
together, these findings demonstrate that H2A.Z functions 
cooperatively with WDR5 and BPTF to stimulate the 
expression of tumorigenic genes. 

H2A.Z is essential for H3K4me and WDR5/BPTF recruitment 

To check whether the cooperative functions of WDR5 
and BPTF in H2A.Z-activated genes reflect their targeted 
recruitment, we next investigated their localization at 
H2A.Z target genes by ChIP assays. Cross-linked chro- 
matin was isolated from control cells and cells depleted 
of H2A.Z, WDR5 or BPTF, and the precipitated DNA was 
amplified by qPCR using primers specific for upstream 
region (UR), transcription start site (TSS) and coding 
region (CR). In agreement with the ChlP-seq data, we 
detected high H2A.Z signals around the TSS, but much 
reduced signals upstream and downstream of the TSS 
in BIRC3 and CCND1 genes (Figure 6A, B, H2A.Z). In 
determining regions occupied by WDR5 and BPTF, we 
found that their binding was near the TSS and coding 
region (WDR5 and BPTF). Consistent with our nucleo- 
some purification studies and previous reports linking 
H3K4me2/me3 to the recruitment of WDR5 and BPTF 
[26-28], the distribution of H3K4me2/me3 at the loci was 
similar to those of WDR5 and BPTF (H3K4me2/me3). 

These data raised the question of whether colocalization 
of WDR5 and BPTF at the target genes is dependent on 
H2A.Z-promoted H3K4me2/me3. In fact, the predicted 
effects of H2A.Z were confirmed by the observation 
that the levels of H3K4me2/me3, WDR5 and BPTF at 
the target loci were reduced in response to H2A.Z 
knockdown. On the contrary, knockdown of WDR5 or 
BPTF did not have substantial effects on H2A.Z occupancy 
around the TSS, placing WDR5 and BPTF in later steps of 
chromatin remodeling. Moreover, H3K4me3 and BPTF 
localization at the target genes decreased sharply upon 
WDR5 knockdown, suggesting that BPTF is recruited 
to the target genes via its physical recognition of WDR5- 
facilitated H3K4me3. Because the conversion of H3K4me2 
to H3K4me3 is attenuated in WDR5-depleted cells, we 
detected a sharp increase in H3K4me2 over the TSS 
and CR of the target genes. As expected from the qRT-PCR 
results (Figure 5A), ChIP analysis showed no enrichments 
of H2A.Z, WDR5 and BPTF at the CASP8 gene, which 
did not alter after individual depletion of H2A.Z, WDR5, 



and BPTF (Additional file 1: Figure S6). In view of the 
importance of H3K4me2/me3 for WDR5/BPTF recruit- 
ment, we also suspected that H2A.Z target genes may 
be selectively enriched for H3K4me2/me3, WDR5 and 
BPTF after ectopic H2A.Z expression in the normal 
cells. Indeed, when stably expressed in the UROtsa cells, 
ectopic H2A.Z produced a distinct increase in the levels of 
H3K4me2/me3 at the target genes and thereby facilitated 
the recruitments of WDR5 and BPTF (Additional file 1: 
Figures S7A and S7B). Consistent with these data, ectopic 
expression of H2A.Z in the normal cells increased target 
gene expression and cell migration (Additional file 1: 
Figures S7C and S7D). These results strongly argue that 
H2A.Z, WDR5 and BPTF function coordinately in regulat- 
ing growth-controlling genes and their concerted actions 
are linked to H3K4me2/me3. 

Discussion 

Despite significant efforts to understand altered expression 
and functional deregulation of H2A.Z in human cancer, 
we are still in the early stages of mechanistic studies for its 
contribution toward cancer development. Given the close 
correlation between H2A.Z exchange and gene activation 
in the context of chromatin, H2A.Z has been implicated 
in upregulating the network of genes during carcinogen- 
esis. To gain deeper understanding of oncogenic effects of 
H2A.Z, we analyzed H2A.Z localization and function in 
bladder cancer cells by ChlP-seq and gene expression 
array. Our analysis of H2A.Z-enriched genes in bladder 
cancer cells revealed a significant overrepresentation 
for tumorigenic genes in general. By comparing tran- 
scriptional profiles, we found that approximately 28% 
of the H2A.Z-enriched genes are activated, whereas 
only approximately 5% of the genes are repressed, in 
the cancer cells. This is significant because H2A.Z is 
known to regulate transcription positively or negatively. 
However, the changes in the expression and localization 
of H2A.Z in bladder cancer cells mainly correlate with 
gene activation, with the most activated genes tending 
to be related to cell growth regulation. Thus, our data 
favor the mechanism in which the elevated expression and 
redistribution of H2A.Z contribute to selective, rather 
than general, disruption of transcription integrity in 
bladder carcinoma cells. The fact that H2A.Z knockdown 
affects the same sets of genes in cancer cells further 
confirms that H2A.Z is a critical driving force for the 
tumorigenic transcriptional output. This study represents 
the first genome-wide analysis of H2A.Z in bladder cancer 
cells, and establishes H2A.Z as a key component whose 
enrichment near the TSS drives expression of major cell 
cycle regulatory genes capable of facilitating cell growth 
and proliferation. Recent studies in human osteosarcoma 
(U20S) cell line showed that H2A.Z is incorporated into 
the p21 promoter region in p53-dependent manner, 
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immunoprecipitated DNA relative to input DNA. 



and establishes a repressive barrier to p53-mediated 
p21 induction [29]. However, we did not observe any 
significant change in p21 transcription after H2A.Z 
knockdown in LD611 cells (data not shown). Furthermore, 
our ChlP-seq data indicate that H2A.Z is incorporated 



into p21 promoter nucleosomes at low levels in both 
LD611 and UROtsa cells. These results imply that 
H2A.Z-catalyzed oncogenic pathway in bladder cancer 
cells is minimally dependent on a negative behavior of 
H2A.Z toward p21 gene. 
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Although the mechanism by which H2A.Z promotes 
chromatin transcription remains unclear, our interaction 
data demonstrated that the presence of H2A.Z in the 
nucleosome increases the levels of H3K4me2 and 
H3K4me3, marks that are associated with gene activation. 
These results are in good agreement with those of 
recent studies showing that H2A.Z colocalizes with 
active histone modifications in transcription regulatory 
regions [7,10,30-32]. It is also conceivable that H2A.Z- 
induced H3K4me2/me3 at the TSS coincides with the 
recruitment of WDR5 and BPTF, which are essential 
components of the MLL H3K4- me thyltransf erase complex 
and the NURD ATP-dependent remodeling complex, 
respectively [33-35]. This suggests that H2A.Z is an 
integral part of histone modification-mediated factor 
recruitment processes and has a direct role in resetting 
transcription competency of growth-regulatory genes in 
bladder cancer. A recent study reported copurification 
of WDR5 with H2A.Z nucleosomes, but BPTF was 
absent in the purified H2A.Z nucleosomes [30]. Clearly, 
these differences could be due to different purification 
approaches (for example, single-step purification using the 
Flag-tag moiety in the previous study versus two-step 
purification using Flag and HA tags in the current study) 
and different wash conditions (150 mM NaCl versus 300 
mM NaCl). We speculate that the weak binding of BPTF 
to H2A nucleosomes is not preserved under our stringent 
purification conditions. It is also noteworthy that H2A.Z 
nucleosomes isolated from bladder cancer cells display the 
levels of the two repressive histone marks, H3K9me3 and 
H3K27me3, similar to those of H2A nucleosomes. Thus, 
although further work will be required to gain a complete 
understanding of how H2A.Z functionally interacts with 
various histone modifications, these findings nevertheless 
illustrate that H2A.Z overexpression in bladder cancer 
cells is likely to act in concert with active histone modi- 
fications which establish constitutive expression of the 
target genes. 

Conclusions 

We show that H2A.Z is overexpressed in bladder cancer 
cells and establishes the active states of oncogenic tran- 
scription. H2A.Z enrichments at target genes increase 
the levels of H3K4me2 and H3K4me3, which in turn 
facilitate the recruitments of WDR5 and BPTF. These 
results unveil a specific contribution of H2A.Z to chro- 
matin remodeling in bladder cancer cells and emphasize 
the importance of accurate levels of H2A.Z expression to 
prevent the onset of tumorigenesis. 

Methods 

Cell culture, constructs and antibodies 

LD611 cells were maintained in Dulbeccos modified 
Eagles medium (DMEM) supplemented with 10% fetal 



bovine serum (FBS). UROtsa cells were cultured in 
DMEM containing 2 mg/ml glucose and 5% FBS. For 
mammalian expression of H2A and H2A.Z, their cDNAs 
were amplified by PCR and ligated into the Nhel and 
BamHI sites of pIRES expression vector containing 5' 
Flag and HA coding sequences. Antibodies used in this 
study are as follows: H2A, H2A.Z, H3K4mel and WDR5 
antibodies from Abeam (Cambridge, UK), Actin, HA and 
Flag antibodies from Sigma- Aldrich (St Louis, MO, USA), 
BPTF antibody from Bethyl Laboratories (Montgomery, 
TX, USA), H3K4me3 and H3K9me3 antibodies from 
Active Motif (Carlsbad, CA, USA), and H3K27me2 and 
H3K27me3 antibodies from Merck Millipore (Darmstadt, 
Germany). 

Tissue microarrays and immunohistochemistry 

The expression patterns of H2A.Z were determined in 
human bladder tissue microarrays, which include 16 
cases of bladder malignant tissues and eight cases of 
normal bladder tissues (US Biomax, Rockville, MA, 
USA). The tissues were deparaffinized, treated with 3% 
hydrogen peroxide to block endogenous peroxide, 
and subjected to antigen retrieval with citrate buffer 
(10 mM, pH 6.0). The tissue sections were blocked in 5% 
normal goat serum in Tris-buffered saline-Triton X-100 
(TBST, 50 mM Tris-HCl, pH 7.5, 150 mM NaCl, and 
0.3% Triton X-100) for 30 min at room temperature 
and incubated with H2A.Z antibody at 4°C overnight. 
Immunodetection was performed using Vectastain 
ABC kit (Vector Laboratories, Burlingame, CA, USA) 
according to manufacturers recommendations. DAB 
(Vector Laboratories) and hematoxylin (Sigma- Aldrich) 
were used for color development and counterstaining. 
The staining intensity was evaluated by semiquantitative 
immunohistochemical assessment, and divided into four 
categories: negative staining (less than 20% of cells), 
weak staining (21 to 50% of cells), moderate staining 
(51 to 80% of cells) and strong staining (more than 
80% of cells). 

ChlP-seq 

H2A.Z ChIP for UROtsa normal and LD611 bladder 
cancer cell lines (1 x 10 8 cells) was carried out following 
a standard protocol as described previously [36]. Briefly, 
cells were fixed for 10 min in 1% formaldehyde, and 
chromatin was sheared with a bioruptor (Diagenode, 
Denville, NJ, USA) to an average size of 200 to 500 bp. 
Lysates were clarified by centrifugation at 20,000 x g for 
15 min at 4°C, and incubated with anti-H2A.Z antibody 
overnight. Protein A/G agarose beads were added to 
protein complexes for 60 min at 4°C. Immunoprecipi- 
tates were washed three times in wash buffer (20 mM 
Tris-HCl, pH 8.0, 300 mM NaCl, 2 mM EDTA, 0.1% 
SDS, and 1% Triton X-100), once in LiCl buffer (10 mM 
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Tris-HCl, pH 8.0, 1% deoxycholate, 1 mM EDTA, 0.25 
M LiCl, and 1% NP-40), twice in TE buffer (10 mM 
Tris-HCl, pH 8.0 and 1 mM EDTA) and then eluted in 
elution buffer (1% SDS and 0.1 M NaHC0 3 ). Crosslinking 
was reversed by overnight incubation at 65°C. After RNase 
A (200 (ig/ml) and proteinase K (200 (ig/ml) digestions, 
ChIP DNA was purified by phenol-chloroform extraction 
and ethanol precipitation. The high-throughput DNA 
sequencing was performed using Illumina HiSeq 2000 
(Illumina, San Diego, CA, USA) at the University of 
Southern California (USC) Genomic Core. 

ChlP-seq data analysis 

Approximately 17 and 13 million sequences were gener- 
ated for the normal and cancer cells, respectively. MAQ 
was used to generate 14 and 11 million aligned reads to 
hgl8 reference genome (mapping quality score >20) [37]. 
One of the duplicated sequences was kept to minimize the 
artifacts of PCR amplification. Since we used single-end 
sequencing, each read was extended in the sequencing 
orientation to a total of 250 bases to infer the coverage 
at each genomic position. The genome was divided into 
non-overlapping windows, and aligned read was consid- 
ered to be within a window of the midpoint of its estimated 
fragment. Mid-points in each window were counted, and 
empirical distributions of windows counts were created. 
The Java-based EpiChIP software was used to identify 
optimal sequence window with respect to gene coordi- 
nates for analysis of H2A.Z enrichment at all Reference 
Sequence (RefSeq) genes. In this analysis, normalized 
ChlP-seq signal for each gene (area below the peak and 
within window) was quantified to yield NLCS (normalized 
locus-specific chromatin state) value by using expectation 
and maximization algorithm and a curve fitting model 
of Baysen information criteria as previously described 
[24,38]. The mapping output files were also converted 
to browser-extensible data (BED) files. For visualization, 
wiggle tracks were generated by computing mean read 
density over 25 bp bins of human genome with aligned 
and filtered reads from UROtsa and LD611 ChlP-seq data. 
Wiggle tracks were visualized in the IGV (Integrated 
Genomic Viewer) [39]. K-means clustering was performed 
on ChlP-seq tag density around 3 kb regions across the 
TSS in R statistical software (http://www.r-project.org). 

Gene expression microarray 

Total RNA was isolated from LD611 and UROtsa cells 
using the TRIzol reagent (Invitrogen, Carlsbad, CA, USA). 
Gene expression microarray experiments were conducted 
using a whole-genome expression array (Human HT-12 
v4 Expression BeadChip; Illumina). This high-density 
oligonucleotide array chip targets more than 47,000 
probe sequences derived from National Center for Bio- 
technology Information Reference Sequence (NCBI) 



RefSeq Release 38 (November 7, 2009) and other sources. 
Data were background-corrected using Genome Studio 
(Illumina) and were quantile-normalized. Genes showing 
detection P <0.05 with fold change >1.5 were selected as 
suggested by the manufacturer (Illumina). Genes which 
are differentially associated with H2A.Z in LD611 cancer 
cells were functionally analyzed in the context of gene 
ontology and molecular networks by using Ingenuity 
pathway software (IPA; www.ingenuity.com). Differen- 
tially regulated genes were categorized into various 
functional groups (threshold P <0.05) and mapped to 
genetic networks. 

RNA interference, qRT-PCR, and ChIP 

For shRNA-based knockdown, DNA oligonucleotides 
encoding shRNAs specific for H2A.Z mRNA (5'-GC 
TTCAAAGAAGCTATTGATT-3 ), WDR5 (5 -GAG AG 
TGGCTGGCAAGTTCAT-3 ') and BPTF mRNA 
(5 ' -CGCCACTAACAGAGAAGGATT-3 ') were annealed 
and ligated into the lentiviral expression vector pLKO.l 
(Addgene, Cambridge, MA, USA). Lentivirus particles 
were generated in 293T cells by co-transfecting plasmids 
encoding VSV-G, NL-BH and the shRNAs. For H2A.Z, 
WDR5 and BPTF knockdown, LD611 cells were infected 
with these viruses and selected with puromycin (2 (ig/ml) 
for two weeks. For qRT-PCR, total RNA (2 \ig) was reverse 
transcribed using the iScript cDNA Synthesis Kit (Bio-Rad 
Laboratories, Hercules, CA, USA) and the PerfeCta™ 
SYBR Green FastMix (Quanta Biosciences Gaithersburg, 
MD, USA) as described [36]. All samples were run in 
triplicate, and results were averaged. Sequences of the 
primers used for qRT-PCR are listed in Additional file 6: 
Table S5. For conventional ChIP assays, LD611 cells were 
harvested, fixed with 1% formaldehyde and processed for 
immunoprecipitation using antibodies against H2A.Z, 
WDR5, BPTF, and H3K4mel/me2/me3 as recently 
described [36]. qPCR was performed in triplicate using 
the primers listed in Additional file 6: Table S5. 

Nucleosome purification and protein identification 

Nuclei were prepared from LD611 cells expressing Flag 
and HA- tagged H2A or H2A.Z and were digested with 
micrococcal nuclease to produce predominantly mono- 
nucleosomes. Mononucleosomes containing ectopic H2A 
or H2A.Z were subjected to anti-Flag M2 agarose. After 
extensive washing with wash buffer (20 mM HEPES-KOH, 
pH 7.9, 0.5 mM EDTA, 300 mM NaCl, 1 mM dithio- 
threitol, 10% glycerol, and 0.1% Nonidet P-40), the eluates 
were further purified using anti-HA agarose affinity chro- 
matography. The imunoprecipitated nucleosome samples 
were analyzed by liquid chromatography-tandem mass 
spectrometry (LC-MS/MS) as described [40]. 
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Cell proliferation and migration assays 

Cell proliferation was assessed by the 3-4, 5- (dimethyl- 
thyazol-2-yl)-2, 5-diphenyltetrazolium (MTT) assay as 
previously described [36]. In brief, cells were seeded in 
24-well tissue culture plates at a density of 2 x 10 4 and 
treated with the MTT labeling reagent (0.5 mg/ml) at 
37°C for 1 h. The blue MTT formazan precipitate was 
dissolved with the MTT solvent (0.2 ml) and measured 
at a wavelength of 570 nm using a microplate reader 
(Bio-Rad). To measure the cell migration activity, Transwell 
migration assays were performed using 5 x 10 4 cells. Cells 
suspended in serum-free medium were plated on the 
Transwell filters (8-um pore size; Corning, Corning, NY, 
USA), and DMEM supplemented with 10% FBS was used 
as chemoattractant. Cells were stained with 0.5% crystal 
violet and quantified after migration for 18 h. 

Additional files 



Additional file 1: Gene dysregulation by histone variant H2A.Z in 
bladder cancer. 
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expression array. 

Additional file 5: Table S4. H2A and H2A.Z nucleosomes interacting 
proteins. 

Additional file 6: Table S5. Primer sequences for RT-PCR and ChlP. 



Competing interests 

The authors declare that they have no competing interests. 
Authors' contributions 

KK and WA conceived and designed the experiments. VP performed 
computational analyses. KK, JC, KH and JK carried out the experiments. 
PL provided reagents and technical assistance. KK, VP and WA wrote the 
manuscript. All authors read and approved the final manuscript. 

Acknowledgements 

We thank Peter Jones for the UROtsa and LD61 1 cell lines. 
Author details 

department of Biochemistry and Molecular Biology, Norris Comprehensive 
Cancer Center, 1450 Biggy Street, Los Angeles, CA 90033, USA. 
2 Bioinformatics Core and Division of Hematology, Norris Comprehensive 
Cancer Center, 1450 Biggy Street, Los Angeles, CA 90033, USA. 3 USC 
Epigenome Center, University of Southern California Keck School of 
Medicine, 1450 Biggy Street, Los Angeles, CA 90033, USA. Research Center, 
Dongnam Institute of Radiological and Medical Sciences, 40 Jwadong-gil, 
Gijang-gun, Busan 619-953, South Korea. 

Received: 6 June 2013 Accepted: 27 September 2013 
Published: 16 October 2013 

References 

1. Luger K, Mader AW, Richmond RK, Sargent DF, Richmond TJ: Crystal 
structure of the nucleosome core particle at 2.8 A resolution. 
Nature 1997, 389:251-260. 

2. Kusch T, Workman JL: Histone variants and complexes involved in their 
exchange. Subcell Biochem 2007, 41:91-109. 

3. Kouzarides T: Chromatin modifications and their function. Cell 2007, 
128:693-705. 



4. Talbert PB, Henikoff S: Histone variants-ancient wrap artists of the 
epigenome. Nat Rev Mol Cell Biol 2010, 1 1 :264-275. 

5. Wang GG, Allis CD, Chi P: Chromatin remodeling and cancer, Part II: 
ATP-dependent chromatin remodeling. Trends Mol Med 2007, 13:373-380. 

6. Zlatanova J, Thakar A: H2A.Z: view from the top. Structure 2008, 1 6:1 66-1 79. 

7. Barski A, Cuddapah S, Cui K, Roh TV, Schones DE, Wang Z, Wei G, Chepelev I, 
Zhao K: High-resolution profiling of histone methylations in the human 
genome. Cell 2007, 1 29:823-837. 

8. Guillemette B, Bataille AR, Gevry N, Adam M, Blanchette M, Robert F, 
Gaudreau L: Variant histone H2A.Z is globally localized to the promoters 
of inactive yeast genes and regulates nucleosome positioning. 

PLoS Biol 2005, 3:e384. 

9. Mavrich TN, Jiang C, loshikhes IP, Li X, Venters BJ, Zanton SJ, Tomsho LP, 
Qi J, Glaser RL, Schuster SC, Gilmour DS, Albert I, Pugh BF: Nucleosome 
organization in the Drosophila genome. Nature 2008, 453:358-362. 

10. Raisner RM, Hartley PD, Meneghini MD, Bao MZ, Liu CL, Schreiber SL, Rando OJ, 
Madhani HD: Histone variant H2A.Z marks the 5' ends of both active and 
inactive genes in euchromatin. Cell 2005, 123:233-248. 

1 1 . Whittle CM, McClinic KN, Ercan S, Zhang X, Green RD, Kelly WG, Lieb JD: 
The genomic distribution and function of histone variant HTZ-1 during 
C. elegans embryogenesis. PLoS Genet 2008, 4:e1000187. 

12. Gallant-Behm CL, Ramsey MR, Bensard CL, Nojek I, Tran J, Liu M, Ellisen LW, 
Espinosa JM: DeltaNp63alpha represses anti-proliferative genes via H2A.Z 
deposition. Genes Dev 2012, 26:2325-2336. 

13. Yang X, Noushmehr H, Han H, Andreu-Vieyra C, Liang G, Jones PA: Gene 
reactivation by 5-aza-2'-deoxycytidine-induced demethylation requires 
SRCAP-mediated H2A.Z insertion to establish nucleosome depleted 
regions. PLoS Genet 2012, 8:e 1002604. 

14. Zhang H, Roberts DN, Cairns BR: Genome-wide dynamics of Htz1, a 
histone H2A variant that poises repressed/basal promoters for activation 
through histone loss. Cell 2005, 123:219-231. 

15. Zilberman D, Coleman-Derr D, Ballinger T, Henikoff S: Histone H2A.Z and 
DNA methylation are mutually antagonistic chromatin marks. 
Nature 2008, 456:125-129. 

16. Henikoff S, Henikoff JG, Sakai A, Loeb GB, Ahmad K: Genome-wide profiling 
of salt fractions maps physical properties of chromatin. Genome Res 2009, 
19:460-469. 

17. Tolstorukov MY, Kharchenko PV, Goldman JA, Kingston RE, Park PJ: 
Comparative analysis of H2A.Z nucleosome organization in the human 
and yeast genomes. Genome Res 2009, 19:967-977. 

18. Suto RK, Clarkson MJ, Tremethick DJ, Luger K: Crystal structure of a 
nucleosome core particle containing the variant histone H2A.Z. 
Nat Struct Biol 2000, 7:1 1 21 -1 1 24. 

19. Weber CM, Henikoff JG, Henikoff S: H2A.Z nucleosomes enriched over 
active genes are homotypic. Nat Struct Mol Biol 2010, 17:1500-1507. 

20. Conerly ML, Teves SS, Diolaiti D, Ulrich M, Eisenman RN, Henikoff S: 
Changes in H2A.Z occupancy and DNA methylation during B-cell 
lymphomagenesis. Genome Res 2010, 20:1383-1390. 

21. Hua S, Kallen CB, Dhar R, Baquero MT, Mason CE, Russell BA, Shah PK, Liu J, 
Khramtsov A, Tretiakova MS, Krausz TN, Olopade Ol, Rimm DL, White KP: 
Genomic analysis of estrogen cascade reveals histone variant H2A.Z 
associated with breast cancer progression. Mol Syst Biol 2008, 4:188. 

22. Svotelis A, Gevry N, Grondin G, Gaudreau L: H2A.Z overexpression promotes 
cellular proliferation of breast cancer cells. Cell Cycle 2010, 9:364-370. 

23. Valdes-Mora F, Song JZ, Statham AL, Strbenac D, Robinson MD, Nair SS, 
Patterson Kl, Tremethick DJ, Stirzaker C, Clark SJ: Acetylation of H2A.Z is a 
key epigenetic modification associated with gene deregulation and 
epigenetic remodeling in cancer. Genome Res 2012, 22:307-321. 

24. Hebenstreit D, Gu M, Haider S, Turner DJ, Lio P, Teichmann SA: EpiChIP: 
gene-by-gene quantification of epigenetic modification levels. 
Nucleic Acids Res 201 1, 39:e27. 

25. Choi J, Heo K, An W: Cooperative action of TIP48 and TIP49 in H2A.Z 
exchange catalyzed by acetylation of nucleosomal H2A. Nucleic Acids Res 
2009, 37:5993-6007. 

26. Li H, llin S, Wang W, Duncan EM, Wysocka J, Allis CD, Patel DJ: Molecular 
basis for site-specific read-out of histone H3K4me3 by the BPTF PHD 
finger of NURF. Nature 2006, 442:91-95. 

27. Wysocka J, Swigut T, Milne TA, Dou Y, Zhang X, Burlingame AL, Roeder RG, 
Brivanlou AH, Allis CD: WDR5 associates with histone H3 methylated at 
K4 and is essential for H3 K4 methylation and vertebrate development. 
Ce// 2005, 121:859-872. 



Kim et a I. Epigenetics & Chromatin 2013, 6:34 
http://www.epigeneticsandchromatin.eom/content/6/1/34 



Page 13 of 13 



28. Wysocka J, Swigut T, Xiao H, Milne TA, Kwon SY, Landry J, Kauer M, Tackett 
AJ, Chait BT, Badenhorst P, Wu C, Allis CD: A PHD finger of NURF couples 
histone H3 lysine 4 trimethylation with chromatin remodelling. 
Nature 2006, 442:86-90. 

29. Gevry N, Chan HM, Laflamme L, Livingston DM, Gaudreau L: 
p21 transcription is regulated by differential localization of 
histone H2A.Z. Genes Dev 2007, 21:1869-1881. 

30. Draker R, Ng MK, Sarcinella E, Ignatchenko V, Kislinger T, Cheung P: 

A combination of H2A.Z and H4 acetylation recruits Brd2 to chromatin 
during transcriptional activation. PLoS Genet 2012, 8:e1 003047. 

31 . Pandey R, Dou Y: H2A.Z sets the stage in ESCs. Cell Stem Cell 201 3, 
12:143-144. 

32. Watanabe S, Radman-Livaja M, Rando OJ, Peterson CL: A histone 
acetylation switch regulates H2A.Z deposition by the SWR-C remodeling 
enzyme. Science 2013, 340:195-199. 

33. Dou Y, Milne TA, Ruthenburg AJ, Lee S, Lee JW, Verdine GL, Allis CD, Roeder RG: 
Regulation of MLL1 H3K4 methyltransferase activity by its core components. 
Nat Struct Mol Biol 2006, 1 3:71 3-71 9. 

34. Dou Y, Milne TA, Tackett AJ, Smith ER, Fukuda A, Wysocka J, Allis CD, Chait BT, 
Hess JL, Roeder RG: Physical association and coordinate function of the H3 
K4 methyltransferase MLL1 and the H4 K16 acetyltransferase MOF. 
Ce//2005, 121:873-885. 

35. Xue Y, Wong J, Moreno GT, Young MK, Cote J, Wang W: NURD, a novel 
complex with both ATP-dependent chromatin-remodeling and histone 
deacetylase activities. Mol Cell 1998, 2:851-861. 

36. Kim K, Heo K, Choi J, Jackson S, Kim H, Xiong Y, An W: Vpr-binding protein 
antagonizes p53-mediated transcription via direct interaction with H3 
tail. Mol Cell Biol 2012, 32:783-796. 

37. Li H, Ruan J, Durbin R: Mapping short DNA sequencing reads and calling 
variants using mapping quality scores. Genome Res 2008, 18:1851-1858. 

38. Hebenstreit D, Teichmann SA: Analysis and simulation of gene expression 
profiles in pure and mixed cell populations. Phys Biol 201 1, 8:035013. 

39. Robinson JT, Thorvaldsdottir H, Winckler W, Guttman M, Lander ES, Getz G, 
Mesirov JP: Integrative genomics viewer. Nat Biotechnol 201 1, 29:24-26. 

40. Heo K, Kim JS, Kim K, Kim H, Choi J, Yang K, An W: Cell-penetrating H4 tail 
peptides potentiate p53-mediated transactivation via inhibition of G9a 
and HDAC1. Oncogene 2013, 32:2510-2520. 



doi:1 0.1 186/1 756-8935-6-34 

Cite this article as: Kim et a I.: Gene dysregulation by histone variant 
H2A.Z in bladder cancer. Epigenetics & Chromatin 2013 6:34. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at f~\ RiftMM i rpntral 

www.biomedcentral.com/submit \^ «««"™a central 



